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We investigate transport properties of one-dimensional fermionic tight binding models featuring 
nearest and next-nearest neighbor hopping, where the fermions are additionally subject to a weak 
short range mutual interaction. To this end we employ a pertinent approach which allows for 
a mapping of the underlying Schrodinger dynamics onto an adequate linear quantum Boltzmann 
equation. This approach is based on a suitable projection operator method. From this Boltzmann 
equation we are able to numerically obtain diffusion coefficients in the case of non-vanishing next- 
nearest neighbor hopping, i.e., the non-integrable case, whereas the diffusion coefficient diverges 
without next-nearest neighbor hopping. For the latter case we analytically investigate the decay 
behavior of the current with the result that arbitrarily small parts of the current relax arbitrarily 
slowly which suggests anomalous diffusive transport behavior within the scope of our approach. 

PACS numbers: 72.10.-d, 05.60.Gg, 05.70.Ln, 75.76.+J 



I. INTRODUCTION 

Considerable effort has been dedicated to the investiga- 
tion of the transport behavior of one-dimensional quan- 
tum wires. While it is accepted that the transport behav- 
ior of the electrons will be normal (diffusive) if the elec- 
trons are coupled to the phonons of the supporting lat- 
tice [l[ , there is a discussion on whether or not electron- 
electron interaction alone will render the transport dif- 
fusive, i.e., non-ballistic. If the electrons are subject to 
a periodic potential, as implied by any standard tight 
binding model the total electron momentum is not nec- 
essarily conserved and one would expect diffusive trans- 
port. However, 1-d spinless fermionic models with near- 
est neighbor hopping and mutual interactions are acces- 
sible by a Bethe ansatz, which means that very many 
local conserved quantities exist, thus they are called "in- 
tegrable" (2J- As follows from the Mazur inequality the 
transport will be ballistic if there is a significant overlap 
of the current with any conserved quantity thus leading 
to a finite Drude weight [3, lJ ■ Such a finite Drude weight 
has been found for zero temperature and small interac- 
tions (gapless regime) for all fillings J5|- A finite Drude 
weight has also been found for finite temperatures and 
all fillings but half filling [3|, |j| . Recently there have been 
papers arguing in both directions for half filling and finite 
temperature (in the gapless regime) : in favor of diffusive 
transport (in spite of integrability) [Q-Q and in favor of 
ballistic transport |9l-[l6j. 

Furthermore, the question how transport behavior 
changes if integrability is continuously broken has been 
addressed. In the low temperature regime an argu- 
ment based on bosonization shows that at least two 
non-commuting umklapp processes are required for non- 
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ballistic transport [ljj. In all temperature regimes dif- 
fusion coefficients may scale as A~ 4 (as opposed to A~ 2 ) 
with A quantifying the integrability breaking term if the 
perturbation allows for the construction of new approxi- 
mately conserved quantities [18| . 

In the paper at hand we investigate a 1-d model of 
spinless fermions with nearest neighbor hopping (nn), 
next-nearest neighbor hopping (mm) and a short ranged, 
perturbatively weak particle-particle interaction. In this 
model the nnn-term is the one that breaks integrabil- 
ity. We address transport behavior by (approximately) 
mapping the quantunr dynamics onto a linear Boltzmann 
equation via a projection operator technique [19 '] . For 
simplicity all calculations are done at infinite temper- 
ature. Within this framework we find concrete trans- 
port coefficients and evidence that they diverge as the 
nnn-term goes to zero. This suggests non-diffusive trans- 
port in the gapless regime at all temperatures. The time 
scales, however, at which this " high- mobility" behavior 
emerges become infinitely long. Results from (approx- 
imate) projection techniques may be wrong if the pro- 
jected subspace does not include all relevant slow hydro- 
dynamic modes. As may be inferred from a paper by 
Belitz [20], however, projective results may at least be 
viewed as reliable lower bounds to the transport coeffi- 
cients. Furthermore, in our projection we do not only 
keep, say, the particle and the energy current but all in- 
dividual occupation numbers of momentum modes. 



II. INTRODUCTION OF THE MODEL 

In this paper we investigate 1-d models of weakly mu- 
tually interacting spinless fermions. I.e., we consider 
fermions on a periodic 1-d crystal lattice. The total sys- 
tem consists of a hopping model, which describes the 
non-interacting fermions, and a nearest neighbor inter- 
action term. The Hamiltonian reads 



H = J 



N 



+Aala n al l+1 a n+1 



(1) 



In the hopping term we incorporate nearest neighbor 
hopping determined by the parameter J, according to a 
standard tight binding model, and also allow for some 
next-nearest neighbor hopping measured by the param- 
eter b. A corresponds to the interaction strength. N is 
the number of lattice sites. One may diagonalize the non- 
interacting system via Fourier transformation and arrives 
at 



* = £ 



t 1 >^ W(q) t t 

e k a\.a k + - }_^ —fi~ a k+ q a i- q a kai , 



k k,l,q 

e k = J(cos{k) + bcos(2k)) , W(q) = -2JA C - iq . (2) 

This Hamiltonian describes an interacting quantum 
gas model with the dispersion relation of the non- 
interacting fermions e k and the interaction parameter 
W{q). (The interaction of course obeys quasi momen- 
tum conservation with the quasi momentum transfer q.) 
For vanishing next-nearest neighbor hopping (b = 0), e k 
corresponds to a cosine band, as usually obtained from a 
tight binding model. 

This model may describe, e.g., weakly interacting elec- 
trons on some 1-d atomic wire, possibly deposited on 
some substrate, or, e.g., as further discussed below, in- 
teracting spin systems, like for example an anisotropic 
spin 1/2 Heisenberg chain (by implementing some map- 
ping of spins onto spinless fermions via Jordan- Wigner- 
transformation) . 



III. TRANSPORT AND DIFFUSION 
COEFFICIENT 

In this section we investigate the transport behavior 
of the above introduced system by directly following a 
scheme described in [21| • Basically, this perturbative ap- 
proach is composed of two steps. Firstly, we map the 
underlying quantum dynamics of a certain set of "occu- 
pation number deviations from equilibrium" in momen- 
tum space onto some pertinent master equation, which 
may in some sense by interpreted as a scattering term of 
a corresponding linear (ized) Boltzmann equation. Fur- 
ther information about the mapping of quantum dynam- 
ics onto Boltzmann equations may be found in 21( and 
references therein. In the approach at hand the number 
of variables of the master equation is determined by the 
(finite) number of lattice sites N. The occupation num- 
ber deviations from thermal equilibrium are described by 
operators 



A J : = (1 - fjWjOj - fjdja] = a]a 3 



~fi 



(3) 



where fj actually corresponds to the equilibrium Fermi 
distribution. 

Note that throughout this paper we focus on the case 
of infinite temperature, i.e., in the following we assume 

/i = l/2- 

As dynamical variables we consider the time de- 
pendent expectation values of these operators, i.e., 
dj(t) :=Tr{Ajp(t)}, where p(t) is the density operator 
which describes the actual state of the system as result- 
ing from the unitary Schrodinger dynamics. 

Our approach is based on a pertinent projection oper- 
ator technique (for details see [131). Within this frame- 
work one has to define a suitable projection operator, 
which maps the current system's state p(t) onto a den- 
sity matrix-like object that only contains the variables 
of interest (here dj(t)) as time dependent quantities. In 
this context the projection is chosen as 



Vp{t) 



n c q 



E 
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dj(t) 



(4) 



The resulting master equation reads 



<ii (t) = ^2 Rik (*)* (*) - $3 Rki ^ dt (*) ' (5) 

k^i k^i 

where the corresponding rates (T = oo, k ^ i) are given 
by 



Rki(t)= I dr~J2 



2 1 



N 2 



(Re(W(* - fc)) - Re(W{k - I))) 2 



cos[-(e, + si - e k - e^k+i )r] 

-i(Re(IF(^-^))-Re(VF(fc-0)) 2 ^ 
coB[\{e k+ e i -e i+h _ l -e l )r]. 



(6) 



I.e., every element of the rate matrix (]6|) is given by an 
integral over a certain correlation function. We especially 
intend to evaluate (|6]) in the limit of small interactions 
and large time t. 

Our approach to transport applies if all rates become 
constant in this limit for times larger than some correla- 
tion times T c jk- If this is the case, we may view the rate 
matrix R^i as the collision term of a Boltzmann equa- 
tion from which a diffusion coefficient may be obtained 
using a type of Chapman- Enskog-like approach, see [21| . 
We omit any details here for brevity. We find for the 
diffusion coefficient (to the order of A" 2 and for T = oo) 



D 



1 
'N 



V k R k iVi 



(7) 



where R l is the inverse matrix to R neglecting the 
eigenspace corresponding to the equilibrium distribution 



fj. Vi is the momentum dependent velocity/current 
vector, the individual velocities are given by the slope 
of the dispersion relation of the non-interacting system 
Vi = dei/di. (Hereby we follow the common idea of iden- 
tifying particle velocities with group velocities.) 

Although we consider finite lattices, the obtained dif- 
fusion coefficients have to be robust against upscaling of 
the system in terms of N — > oo in order that our results 
are reasonable for large ("realistic") crystals. 



IV. NUMERICAL RESULTS FOR THE 
DIFFUSION COEFFICIENT FOR b / 

For 6 ^ the analytical evaluation of the long time 
limit of the rates © is not directly accessible. Instead we 
determine Rki(t) numerically for finite time t. (Numer- 
ics, which we omitted here, suggest that the rate matrix 
elements Rki(t) may be expected to become time inde- 
pendent after some correlation times T Ct ik-) Based on 
this time dependent rate matrix we may determine an 
initially also time dependent diffusion coefficient D(t) by 
directly evaluating formula (|7J|, which is our first main 
numerical result. 



9 



We essentially find the following two quantitative re- 
sults. The long time diffusion coefficients become larger 
for smaller 6. This finding clearly demonstrates the quan- 
tum nature of the model and the approach, because 
it seems counter-intuitive in the sense that it does not 
fit the (classical) perspective that an increased hopping 
would lead to a higher diffusion. 




FIG. 2: Long time diffusion coefficients versus next-nearest 
neighbor hopping b. The numerics (points) are weli approx- 
imated by the function 4/V6 (solid line). D(b) diverges for 
6 = 0. 
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FIG. 1: Time dependent diffusion coefficients for system with 
next-nearest neighbor hopping b — 0.2 (white triangles), b = 
0.3 (black triangles), b = 0.4 (white squares), b — 0.5 (black 
squares), b — 0.6 (white circles), b = 0.8 (black circles). N — 
500, A = 0.01. 



Fig. [T] demonstrates these time dependent diffusion co- 
efficients for different next-nearest neighbor hoppings b. 
One finds that for all analyzed b the diffusion coefficients 
D(t) reach plateau levels and therefore become approx- 
imately time independent for large enough times. The 
constancy of D of course indicates that the elements of 
the rate matrix R^i have converged, too. I.e., within this 
framework we are able to numerically determine finite 
long time diffusion coefficients for non-vanishing next- 
nearest neighbor hopping 6^0. 



Additionally, Fig. [2] shows the dependence of the long 
time diffusion coefficients on the next-nearest neighbor 
hopping b. The plot suggests a dependence D{b) oc 6 -1 / 2 
(, although we do not have an analytical proof for this 
type of relation), which would mean that the diffusion 
coefficient diverges for 6 = 0. 

Secondly, the time, after which the diffusion coefficient 
is approximately constant, becomes larger for smaller b. 
I.e., to analyze the interesting limit b — > (see below) 
we have to evaluate longer and longer times to obtain 
the plateau level (, which goes along with larger chain 
lengths N; recall that the results should be robust against 
upscaling of N). 

Since for 6 = the plateau is reached only at infinite 
times and the plateau height is expected to be infinite, 
which is in some sense equivalent, we cannot numerically 
extract finite diffusion coefficients in this case. The find- 
ings may rather suggest that the transport is possibly 
non-diffusive for 6 = 0. Some analytical treatment of 
this case can be found below in Sec. [V] 

At this point this result should be compared to a re- 
cent, related result by Steinigeweg et al. (7J. They con- 
sider the same model class but only for 6 = and employ 
a related projection technique. However, other than in 
the work at hand, they do not keep all individual mo- 
mentum occupation numbers but project onto a single 
variable, namely the current j — X)j u i a / a i- Doing so 
they find a finite diffusion coefficient emerging at short 
time scales for all A including arbitrarily small ones. We 
checked that this scenario gradually transforms into the 



picture we find as more and more observables in addition 
to the current are kept in the projection, i.e., time scales 
become longer and diffusion constants increase. This, 
however, is, according to Belitz 20], not a contradiction 
since results from projections onto subsets of observables 
yield reliable lower bounds to the true transport coef- 
ficients. Thus we conclude, in spite of the findings by 
Steinigeweg et al. that there is no regular diffusive trans- 
port for 6 = 0. 

Concerning the short time/strong interaction behavior 
of the diffusion coefficients we find that the graphs for 
different 6 in Fig. Q] intersect at some time t. It may be 
thinkable that the behavior of the diffusion coefficients 
reverses for small times in the way that the diffusion co- 
efficient for small times becomes larger for larger 6, which 
would be in agreement with (or at least not contradictory 
to) the above mentioned intuitive expectation. However, 
since in our approach the rate matrix is not converged 
out at those times, we cannot assure that the short time 
values are quantitatively correct, they rather serve as an 
indicator for the overall qualitative behavior. For this 
regime also see |22j |. 



V. ANALYTICAL CONSIDERATIONS ON 
TRANSPORT FOR b = 

The special case of vanishing next-nearest neighbor 
hopping b = is of particular interest because the re- 
garded model of spinless fermions ([2J is effectively equiv- 
alent to an anisotropic spin-1/2 Heisenberg chain. The 
corresponding spin Hamiltonian 



corresponding integral, i.e. 



H = J 



N 

E 
=i 



l¥ 
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h.c] +As z n s 



n+l 



(8) 



may be transformed into the Hamiltonian (|TJ) with 
b = via a suitable Jordan- Wigner-transformation. Note 
that we consider here the case of small A, which rather 
corresponds to the case of strong anisotropy in the spin 
picture. 

Recall that in our numerical evaluation in Sec. IIVI we 
have found diverging long time diffusion coefficients for 
6 = which we suggested to interpret as an indicator 
for possibly non-diffusive transport. If this is actually 
the case, one would expect to find a number of conserved 
quantities, at least some of which should have a signifi- 
cant overlap with the current. (Since the model is inte- 
grable if (and only if ) b = (see Sec. HJ), there are a large 
number of conserved quantities for 6 = 0. The transi- 
tion from non-integrability to integrability may here be 
related to the transition from diffusive to non-diffusive 
transport.) 

For 6 = one may analytically determine the rates 
Rki in the limit t — > oo. Assuming that the momentum 
modes are lying densely in /c-space (i.e., assuming long 
chains,) one may substitute the sum over I in © by a 



^ N 



(9) 



and regard the energy terms e+(Z) := (e, + £\ — Ek — 
Si-k+i)) in the first and e_(Z) := (e k +£i -e i+k -i -£j) in 
the second part of the sum as continuous functions of I. 
Routinely exploiting the properties of the sine-function 
one may carry out the time integration in (|6]) and finds 
that in the limit t — >■ oo only the exactly energy con- 
serving scattering processes contribute, i.e., only addends 
featuring (e+ = 0) in the first and (e_ = 0) in the second 
term. 

The second "minus" term yields diverging 5-likc con- 
tributions on lines, with respect to a graphical illustra- 
tion of the scattering operator, defined by j + k = tt and 
j + k = —ir, but only exactly on those lines. For all other 
j, k, even in the direct vicinity of those lines, the "mi- 
nus" term vanishes, i.e., the corresponding contribution 
to the rate matrix from the second term may be written 
as R~ (j, k) oc —5{j + k — n) — S(j + k — n) + S(j — k) for 
t — > oo. This feature may be illustrated by a finite time 
sketch of the total rate matrix Rjk (without diagonal el- 
ements) (see Fig. [3]). One finds large negative matrix 
elements in the vicinity of the above mentioned lines, 
for larger t the negative contributions move closer to the 
lines, which eventually leads to the <5-like structure. 
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FIG. 3: Finite time sketch of the rate matrix Rjk for 6 = 
without diagonal elements, the labels denote the correspond- 
ing momentum mode numbers. There are large (absolute) 
negative scattering matrix elements (here black) between 
modes near lines with j + k = n and j + k = — it, which 
originate from the "minus" term in [6] Parameters: iV = 200, 
J = -1, A = 0.01, T = oo. 

Although these terms are problematic for the determi- 
nation of the complete rate matrix, they yield no contri- 
bution to the diffusion coefficient ([7]) because of symme- 
try arguments (see below) and may therefore be omitted 
in the transport investigations addressed in this paper. 

The complete rate matrix for 6 = features mirror 
symmetry to the band middles, i.e., lines with k = tt/2 
and k = —tt/2, because the cos-function simply passes 
into its negative under inversion on those lines. That 
is, every eigenvector of the rate matrix has to transform 



into itself or its negative under mirror imaging on those 
lines and the eigenvectors may accordingly be classified 
by, say + and — . The matrix therefore separates into 
these symmetry subspaces, i.e., if one chooses a matrix 
representation where all basis vectors feature either + or 
— symmetry, there is no coupling between the two sub- 
spaces, the corresponding matrix elements have to van- 
ish. Note that the current features + symmetry. For any 
+ vector d(j) one finds that d(j) = d(ir — j) for j > and 
d(j) — d{— 7T — j) for j < 0. That is, the multiplication 
of the "minus" part R~ (j, k) with a + vector 



0.02 



d(j) = J dkR'(j,k)d(k) 

= -d(-j-7r)-d(-j + ir)+d(j) (10) 

is equal to 0, since two terms cancel each other, depend- 
ing on whether j is positive or negative, and the remain- 
ing term does not contribute to the integration because 
the corresponding fc-argument is not in the first Brillouin 
zone. Therefore R~ (j, k) gives no contribution to the 
sum evaluated in (|7|) , which means that R~ (j, k) can be 
neglected. 

From the remaining first "plus" term in (|6|) one finally 
obtains from collecting and evaluating all the zeros of the 
respective argument in © 



R + (j,k)<x 



1 (cos(fc - j) + cos(fc + j)) 2 



N 



|sin(fc) -sin(j)| 



(11) 



for t = oo. One recognizes that the rate matrix R + (j, k) 
is here a markedly complicated object, mainly because of 
two reasons. 

There are diverging elements of R + (j, k) in the vicin- 
ity of lines defined by j + k = tx and j + k = —n. For 
larger and larger systems (N — > oo), when the momen- 
tum modes become denser and denser in fc-space, there 
are more and more elements in the direct vicinity of those 
lines, which therefore become arbitrarily large. From this 
one may expect that at least some eigenvalues of R + (j, k) 
diverge in the limit N — !> oo. 

Secondly, there are whole lines in the matrix (|lip de- 
fined by j, k = — 7r/2 and j, k = n/2 where all ma- 
trix elements are 0, indicating conserved quantities. I.e., 
d(j) oc S(j — 7r/2) + S(j + it/2) is a conserved quantity, 
which features a non-singular but finite overlap with the 
current. This causes an infinitely small portion of the 
current to decay infinitely slowly (see Fig|4j) In Q it is 
pointed out that a non-diffusive behavior of the system 
would require a non-local conservation law. Note, how- 
ever, that the above conserved quantity features such a 
spatially non-local structure. 

To take a closer look at the concrete transport behavior 
we evaluate the overlap of the current/velocity vector 
with the eigenvectors of R + (j,k) in dependence of the 
corresponding eigenvalues. 

Firstly, Fig0] shows that the current does not decay 
monoexponentially, although an average relaxation rate 
may be identified. In particular there is a finite overlap 
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FIG. 4: Overlap P n :— (j ■ r n ) 2 of the current j with the 
eigenvectors f„ of the matrix R + (j, k) versus corresponding 
"rescaled" eigenvalue (r n A 2 are eigenvalues of R + (j,k)) for 
time t — > oo. There is a finite overlap with (approximately) 
conserved quantities for r„ « 0, i.e., parts of the current decay 
arbitrarily slowly. Parameters: iV = 800, J = —1, T — oo. 



between the current and conserved (arbitrarily slowly re- 
laxing) quantities, i.e., an infinitely small portion of the 
current does not decay at all. This can explain why we 
are not able to find converging diffusion coefficients ac- 
cording to formula ([7|) and why the numerically calcu- 
lated curve D(t) (cf. Fig. [TJ does not exhibit a plateau 
for 6 = 0. Fig[J]also demonstrates that there are many 
quantities which have no overlap with the current. 

Within the scope of the master equation dynamics §5§ 
the current autocorrelation function C(t) := (j(t)j) de- 
cays multiexponentially, i.e., 



c(t)=j2 p " 



>r n A 2 t 



(12) 



where r„A 2 are eigenvalues of the rate matrix R + (j, k) 
(r n < 0) and P n corresponds to the overlap of the current 
vector with the corresponding eigenvector f n of R + (j, k) 
defined by P n := (j ■ f n ) 2 . For large systems one may 
substitute the sum by an integral 



C(t) = / P(r)e 



rA't 



dr 



(13) 



with continuous eigenvalues rA 2 and the function P(r) 
essentially given by Fig. 0] The long time behavior is 
mainly determined by the regime of small r. Therefore we 
approximate P(r) by the corresponding Taylor expansion 
around r « 



P(r)«P(0) + X;, 



(14) 



with some coefficients c„. Inserting (|14p into (fT3|) one 
finds that for times t larger than 1/A 2 the behavior of 
C(t) is dominated by the term C(t) sa P(0)/(A 2 t) which 



results from the zeroth order term of ((14)) . All other 
contributions decay with t~ 2 or higher orders in t or ex- 
ponentially with t. Although C(t) is essentially a sum of 
exponential functions, it does not behave exponentially 
in this limit. A time dependent diffusion coefficient D(t), 
which may be decomposed as 



D(t) 



l/A 2 f t 

C(t')dt' + / C{t')dt' , 

[i ./l/A 2 



(15) 



is dominated by the second term for times larger than 
l/A 2 . In this limit D(t) behaves logarithmic as 



D(t) oc (l/A 2 )ln(A 2 i) 



(16) 



and consequently does not become constant for t — > oo. 

Also note that the integral over the current- 
autocorrelation function as, e.g., given in (|15|) is directly 
related to the time dependent growth of the mean square 
displacement (MSD) of some initial distribution, the cor- 
responding connection has been established in |23j |. 

Within the scope of our approach the result (fTB"]) would 
lead to the conclusion that the transport is not regu- 
larly diffusive and neither really ballistic (in this case 
one would find D(t) oc t), but rather something in be- 
tween, which may, e.g., be termed "anomalous diffusion" 
(cf. also [24]). This characteristic behavior arises from 
the fact that the function P(r) (cf. Fig. QJ exhibits a 
finite value at r = 0. Ballistic transport would require a 
peak at r = 0. From (fl6|) one may suspect that the MSD 
would grow proportional to tln(t). 



However, our approach is restricted to diffusion coef- 
ficients of the order of A 2 , i.e., there could possibly be 
diffusive behavior scaling with A or higher orders. 



VI. CONCLUSION 

We demonstrated that the quantum dynamics of mo- 
mentum mode occupation numbers in a 1-d quantum 
model for weakly interacting fermions may be mapped 
onto a master equation, the corresponding rate matrix 
can serve as collision term of a linear (ized) Boltzmann 
equation. We numerically found that adequately calcu- 
lated diffusion coefficients become constant in the limit 
of large times and large systems for finite next-nearest 
neighbor hopping. These long time diffusion coefficients 
are larger for smaller next-nearest neighbor hopping. For 
vanishing next-nearest neighbor hopping the diffusion co- 
efficients diverge. In this case we have analytically shown 
that arbitrarily small parts of the current do not re- 
lax because there is a finite overlap between the current 
and (approximately) conserved quantities. This result 
suggests anomalous diffusive transport behavior for the 
model featuring no next-nearest neighbor hopping. 
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